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Abstract 

We deduce Levinson 's theorem in non-relativistic quantum mechanics in one di- 
mension as a sum rule for the spectral density constructed from asymptotic data. 
We assume a self-adjoint hamiltonian which guarantees completeness; the potential 
needs not to be isotropic and a zero-energy resonance is automatically taken into ac- 
count. Peculiarities of this one-dimension case are explained because of the "critical" 
character of the free case u(x) = 0, in the sense that any atractive potential forms 
at least a bound state. We believe this method is more general and direct than the 
usual one in which one proves the theorem first for single wave modes and performs 
analytical continuation. 
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1 Introduction. 



Two generic results in potential scattering stand on their own, and hold with wide gener- 
ality. The first, the optical theorem stems from the fact that the scattered "matter" is 
taken away from the incoming wave, and hence the scattering center "casts a shadow" in 
the forward direction as to produce negative interference with the incoming beam; therefore 
a relation must exist between the total scattering cross section and the forward scattering 
amplitude. Originally proven in partial waves for 3D scattering, the theorem holds with 
much more generality; a simple but very general proof is offered in The theorem can be 
seen as a result of the "completeness relation" in ordinary space (also called orthogonality), 
at a given energy. 

The second result is Levinson 's theorem, which in a way can be seen also as a con- 
sequence of completeness in momentum space. In its primitive form of 1949 Levinson 's 
theorem reads |3j 



for the number rig of bound states of angular momentum £ in a generic central potential 
u(r) = u(\r\) which produces a phase shift Se(k) for scattering with energy E = k 2 ; the 
proof came up as a byproduct of studies on uniqueness of potentials with a given phase shift 
|E]. Formula (JTJ) was rediscovered in 1956 jHj- In 1957 Jauch [6 established the theorem as 
a consequence of the completeness relation for scattering states and set it up in the general 
frame of operator theory. 

The philosophy was that, if a potential generates bound states, there should be a 
relation between them and scattering states as the completeness relation has to be "shared" 
among them. Indeed, Levinson 's theorem is the only relation between bound states and 
scattering states as it can be deduced from inverse scattering theory. Since the work of 
Jauch, many studies followed and we mention in particular the elementary deduction by 
Wellner [7j for s- waves in 3D, which we shall generalize in this paper and the later studies 
of Newton relating the theorem to the inverse problem in 1, 2 or 3 dimensions [Hj, [0] • 

It is considerably more difficult to prove Levinson 's theorem than the optical theorem, 
although both share complementary physical foundations. In this paper we prove the 
theorem as a sum rule for the spectral density, which we take as more fundamental entity; 
we shall work in one dimension with local potential where all the features of the problem 
already show up. Indications for D > 1 will be given at the end of the paper. 

In our one-dimensional problem we shall not assume parity invariance (corresponding to 
central forces in D = 3) nor we shall exclude a zero-energy resonance. We were motivated 
by the spectral density considerations of Niemi and Semenoff ^H] for fermions in solitonic 
backgrounds whereas the wronskian-like technique is adopted from Wellner as stated. The 
plan of the paper is as follows: in sect. 2 we set up the scattering problem in D = 1, 
mainly to motivate notation for direct (left to right) and right to left ("zurdo") scattering, 
unitarity of the S-matrix, etc. In the third section we introduce the spectral density which 
needs to be regularized in a box, but in the definition of relative spectral density the space 
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cutoff can (and will) be removed. As stated, Levinson's theorem will appear clearly as a 
sum rule for the relative spectral density. 

In section fourth we shall handle two simple examples which can be worked out directly, 
namely the solitonic Poschl- Teller potential u(x) = —2 sech 2 x and the delta potential 
u(x) = g S(x); the difference of generic vs. critical potential will be cleared up, as well as 
the one-half factor already noted by Barton [TT], and present in the one-dimensional case 
(another factor should appear in two dimensions) [12] . The next section will exhibit our 
general treatment of the relative spectral density for an arbitrary local potential u = u(x) 
(of course, decaying at x — > ±00 fast enough to allow for scattering). We shall emphasize 
that the spectral density is a "hard datum" (i.e. spectrum-dependent as opposed to the 
rest of potential parameters or "soft data"); the spectral density is given in terms of the 
derivative of the phase of the forward amplitude, this amplitude itself being a hard datum 
as well, explaining a result used routinely in the KdV-like evolution equations. In section 
6 we carry out the momentum integration of the density to produce the general form of 
Levinson's theorem; we comment briefly on the relation with the determinant of the S- 
matrix which suggests an interpretation of the theorem as an index theorem for the bound 
states. In sect. 7 we set up the procedure for an arbitrary dimension, thus generalizing 
the method of Wellner; noncentral and critical potentials (i.e., producing a zero-energy 
resonance) are easily included; finally we make some comments on non-local potentials 
and add some concluding remarks. 



2 Scattering in one-dimension. 

Let 



ip"(x) + k 2 i)(x) = u(x)ip(x) (2) 

be the Schrodinger equation in one dimension for a local (hence real because hermitean) 
potential allowing scattering, i.e. satisfying [T3] 

(1 + x 2 )\u(x)\ dx < 00 (3) 

for positive E = k 2 . The D — 1 scattering has card S° = two modes, direct (incoming 
wave towards the right), with the asymptotic solution: 

< J exp(ikx) + b{k) exp(— ikx) for 1 C , . 

*\ t(k)exp(ikx) for>0 W 

and zurdo scattering: the incoming wave travels towards the left 

r, , J exp(— ikx) + b(k) exp(ikx) for x 3> , , 

V{ x ) > I t(k) exp(-ikx) for x < ' ' 

The amplitudes f(k) = t(k) — 1 and b(k) give the scattering coefficients by the relations 
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a^ = \f(k)\ 2 , ^ = \b(k)\ 2 (6) 

together with 

\t(k)\ 2 + \b(k)\ 2 = l or \f(k)\ 2 + \b(k)\ 2 = a tot = -2Mf(k), (7) 
the optical theorem in one dimension. The S-matrix has two channels only 

*<*>=($! ID < 8 > 

Unitarity S^S = SS^ leads to the important relations ( with a = |a|exp(i$ a ) for any 
amplitude) 

\t(k)\ 2 + \b(k)\ 2 = \t(k)\ 2 + \b(k)\ 2 = \t(k)\ 2 + \Hk)\ 2 = 1 (9) 

and 

t(k)b(k)* + b(k)t(k)* = 0, $[t(k)] + $[t(k)} = 7T + $[&(*;)] + $[6(ife)] (10) 

Now the potential being local (and hermitean of course) is real, so time reversal holds; 
it follows that 

t(k) = t(k) (11) 

whereas an even potential would imply similarly b(k) = b(k), which we do not assume. To 
see (jUJ), notice t{k) is defined as the transition from ki n to k out . Time reversal changes k 
to —k and in to cra£; hence t(k) goes to t(/c). Also we shall take k > for direct scattering 
and k < for zurdo scattering; indeed then t(/c) = —t(—k) for the above reasons. 

We remind that we also use the terms hard data for spectral data, namely the spectral 
density, and soft data for orientation data, to wit, norming constants for bound states and 
the phase of b(k). 



3 The Spectral Density. 

We recall first elementary properties of matrices. Completeness for a diagonalizable finite 
matrix M means 

M = J2™P m or 1| =J2 P m (12) 

m m 

for eigenvalues m and projectors P m ; the second relation is called resolution of the identity. 
In his work on Quantum Mechanics, von Neumann [TJ] extended the classical work of 
Hilbert on integral equations for hermitean unbound operators with continuous spectrum: 
he called hypermaximal (today "self- adjoint") those hermitean operators H which still 
support a resolution of the identity. The resolution reads 



4 



N 



■" poo 

i|« = E p i+/ dp 0") ( 13 ) 

where W is the Hilbert space of states, 1| = l\n is the unit operator, Pj projects to 
the finite or infinite number (N = 0,1,...) of bound states and the continuum integral, 
supposed by simplicity extended from to oo, as in a standard potential problem, means 
projection-valued measures. Notice the unit operator l\ n is bounded but not in the trace 
class. 

The simplest approach to the Levinson's theorem is to state the same resolution for 
the free system H Q 

POO 

l\ n = dpW(Ai) (14) 
Jo 

Substracting in (fT^jl and then taking traces we get, supposing short-range potentials 
which support at most finite number N of bound states 

-JV= r Tr d[PU) - P (o) (a0] (15) 
Jo 

which is, really, the most general (but rather useless) form of the theorem. The idea is now 
to trade the projectors for scattering amplitudes (or phase shifts). We know of course the 
normalized free continuum wavefunctions 

' l l )< k\ x ) = — exp(ikx). (16) 

v27T 



Let us normalize the continuum wavefunctions to 1, instead to V2tt- Then (JT5j) for our 
case really means 

/oo roo 
dk / dx [\M X )\ 2 ~ 1] (17) 
-oo J —oo 

as dP/dk ~ \k >< k\, taking the traces in x-spaces, defining ipk{%) ■=< x\k > and using 
the double degeneracy E = (±k) 2 to extend the integral to ±oo. Expression (fTTj) will be 
our fundamental formula; the /c-integration will be clarified below. 

The main concern is to express (fTTj) in terms of the asymptotic data (phase shifts). 
Define the relative spectral density for the problem as 

/oo 
dx [|^.(x)| 2 -l] (18) 
-oo 

Notice the individual spectral densities diverge, i.e. p/ ree (^) = J 1 dx — oo; only 
relative densities make sense. The idea of the proof is to relate the spectral densities to 
scattering data; as we know that the integral of the spectral density gives Levinson's, 
which express the number N of bound states in terms of the range of the phase shift ( 
— > oo), (see (JTJ), we expect the spectral densities to be given in terms of derivatives of 
the phase shifts; we shall see that this is so. Notice also the relative spectral density is 
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a measure, that is, something under an integral symbol; hence some apparent pathologies 
like delta-type behaviour or box-normalization and removal are perfectly legitimate, and 
not only "heuristic" as some authors state. 



4 Two simple examples: Poschl-Teller Soliton and Delta poten- 
tials. 

As a warming-up exercise, let us compute the spectral densities in two simple cases in 
which the exact continuum wavefunctions are known. When considering the Poschl-Teller 
potential (the standard solitonic potential) given by 

u(x) = ^— , (19) 

cosh x 

there is a single bona fide bound state with energy E — — 1 and continuum states with no 
reflection since the potential is transparent. Indeed, the whole E > wavefunction if)k(x) 
is obtained from the u(x) = case by Darboux's method |15] . If D = d/dx 

ipk( x ) ~ (D — tanh x) exp(ikx) (20) 
With the correct normalization included (so ipk{ x <C 0) ~ exp(ikx)), 

, . , (ik — tanh x) . , . , . 

Mx) = rr— exp zb 21 

ik + 1 

The relative spectral density is therefore 

POO 2 

Prd {k) = / [|^ fc (x)| 2 - 1] dx = — — (22) 



k 2 + l 

The sum rule or fc-integration gives of course 

/oo 
p rel (k) dk = -2tt (23) 
-oo 

so that iV = 1 as expected. There is more to say in this case, e.g. we find a zero-energy 
resonance or "half-bound" state which corresponds to the k = limit, i.e. 

ip k=0 = — tanh x (24) 
As regards the delta potential we have 

u(x) = g 5{x) (25) 

where in principle we leave the sign of g open. As the support of the potential is a point, 
{0}, the solution for x ^ is always asymptotic; there is also no odd wave. For the spectral 
density we compute in this case 
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so that 



/U fCO 
l\Mx)\ 2 -±}dx+ / [\i, k {x)\ 2 -l]dx (26) 
-oo JO 

roc /'OO 

/ [|^(*)| 2 -l]ds= / [|t(A;)| 2 - l}dx = -\b(k)\ 2 L (27) 
Jo Jo 

for L — > oo. This divergence is in fact spurious and cancelled with the x < contribution, 

i.e. 

/° [|?M^)| 2 - l]da; = f [1 + |6(A;)| 2 + 2K{6(A;) exp(-2iib;)} - 1] cfe (28) 

J -co J —CO 



-co 

Now we define 



2 /° $l{b(k) exp(-2ikx)} dx = A + B (29) 

J -co 

so that 

A = 2 f° ^{b(k)} cos 2kx dx (30) 

If we bear in mind that 

/0 rco rco 

cos 2/cx dx — cos 2/cx dx — exp(2ikx) dx = 2jr5(2k) = rr5(k) (31) 
-oo J— CO J— CO 

then we have that 

A = $l{b(0)}n5(k) (32) 

It is the case that 

6(0) = -1 (33) 
for a generic potential, including the delta, because ip k=0 (x) = so 

exp(ikx) + b(k) exp(-ikx) = as k — > (34) 

and therefore 6(0) = —1. The exception (critical potentials) occurs for a zero-energy 
resonance, see later, when 6(0) = 0. On the other hand (L — > +oo) 

B „ r° . , „, N , l-cos2£;L 



3{6(A;)} 

so that 



2 / sin(-2A;a;) cfe = C ° S fc (35) 



B = MMI (36) 
k 

7 



as the oscillatory part cos2kL gives no contribution as a measure when L — > oo. Now for 
the delta potential itself we have (no odd wave) 

/(*) = t(k) - 1 = b(k) = (37) 
thus confirming 6(0) = — 1. So the relative spectral density as a whole is 

Prei(k) = -<*5{k) + — ^ (38) 

Notice the delta piece, which will persist for any generic potential. Also, the dependence 
p re i{k) oc 1/k 2 for k ^> is general as the phase shift itself will fall with 1/k (validity of 
the Born approximation) and we expect p re i{k) oc 5'(k). Now a ^-integration would yield 
Levinson's theorem; taking care to isolate the sign(g) piece we get 

- 2ixN = -7T + sign(g) n (39) 

In other words 

N = 1 = Si f {9) (40) 
which is obviously correct: N = 1(0) for g < (g > 0). 

5 General calculation. 

Now we carry out the general calculation. Starting from 

ip"(x) + k 2 tp(x) = u(x)ip(x) (41) 
we differentiate with respect to k (represented by the dot symbol), i.e. jjj, [IE] 

ip"(x) + 2ki) + k 2 i){x) = u{x)ip{x) (42) 

Why this unusual derivative? Because we expect the spectral density to depend on 
fc-derivatives of the scattering amplitudes (or phases shift) as said before. Next we take 
real and imaginary parts of the wavefunction 

Mx) = ®{Mz)} + S{Mx)} ( 43 ) 
since, the Schrodinger operator being real, each works separately. Now we construct the 
wronskian for, first, 3l{ipk{x)} '■= which satisfies 

&'{x) + 2k<S>(x) + k 2 $(x) = u{x)$(x) (44) 

If we multiply and substract in the usual way (e.g. to get the current) it is the case 
that 
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(x) - $(x)&(x)]' = 2M> 2 (x) (45) 

or (L — > oo eventually) 

sft{J fc } : = - $(x)$'(x)| L L = 2A; ^ $ 2 (x) dx (46) 

which, together with the imaginary contribution, is the crucial result since it allows us to 
express the spectral density in terms of the asymptotic data. Next we define 

A := SR{4} at L, B := %{I k } at L (47) 

C := U{I k } at — L, D := %{I k } at - L (48) 

and 

t(k) = \t(k)\exp(iip t ), b(k) = \b(k)\ exp(iip r ) (49) 
with A := ip t + kL. In doing so 

A = \t(k)\ 2 [k{ip t + L)+ cos A sin A] (50) 

B = \t(k)\ 2 [k((p t + L) - cos A sin A] (51) 
So the total forward contribution is 

A + B = 2k \t(k)\ 2 (<p t + L) (52) 

This is very nice: the factor 2k of ()45|) appears, as well as the derivative of the forward 
phase, (a hard datum; see below) while the L divergence will be spurious. 

The calculation of the backward part is more involved as the wavefunction is 

tj)(x <C 0) = exp(ikx) + b{k) exp(-ikx) (53) 

By repeating the former method, now for the total backward contribution, we get (after 
some cancellations between C and D) 

C + D = AkL — 2\b(k)\ 2 k(<p r - L) + 2b(k) sin(<^ r + 2kL) (54) 
As regards the regularized spectral density pi{k) we have finally A + B — C — D or 

2kp L (k) = 4kL + 2k(p t + 2k\b{k)^{Cp r - <p t ) + 2b(k) sin(^ r + 2kL) (55) 
So the final result is 

- 2tcN = [[ lim (p L (k) - pf{k)} dk (56) 
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where as expected 



pf(k) = f L dx = 2L (57) 

To sum up 

- 2ttN = r [ip t + \b(k) | 2 (0 r - <pt) + h -^r- sin(^ + 2kL)] dk. (58) 

J —CO rC 

As a matter of fact the first term of the integrand would be the relative spectral density. 
What about the back phase and the oscillatory third term? We find that, first, 

r shi^ + 2kL) dk = vr6(0) (59) 

J -co K 

as shown in the Appendix. So the sin(y9 r + 2k L) / k integral is really irS(k) as a distribution. 
This is completely rigorous because we are talking of measures (project ion- valued measures) 
and the delta is itself a measure (not so the delta prime). 

As regards the second term in (J58|) we write the zurdo contribution with k > and 
integrate then from to oo, i.e. 

\b(k)\ 2 (<p r - cp t ) + \b(k)\ 2 (^ r - cp t ) = \b(k)\ 2 (<p r + (p r - 2<p t ) = (60) 

because (see © and (fTUj) ) \b(k)\ = \b(k)\ and 2ip t — ip r — (p r = ir. So the final expression 
for the relative spectral density is 

Prel {k) = if t + nb(0)5{k), k>0 (61) 

see e.g. ^H], ^7]. The first term, derivative of the forward phase, comes by no surprise and 
represents the germ of Levinson's theorem. The second one gives a universal contribution 
since 

b(0) = -1 (generic); 6(0) = (critical) (62) 

as we already discussed. We remark here how the spectral density is a hard function and 
therefore the forward phase, but not the backward one, is a hard datum. Indeed, at least 
for finite-range potentials the forward amplitude can be expressed in terms of the bound 
states plus an integral over the modulus of the reflected amplitudes [TH]. Again the usual 
proof is based on analytic continuation, whereas ours stems directly from the definition of 
hard data as spectral data. 

We can easily compute the value of p re i{k) for large k. From the Born approximation 

1 r°° — % 

f(k) = t(k) - 1 ~ — / u(x) dx = — <u> (63) 
lik J -co Ik 



Hence 



tan^ t (fc) ~ — ~ <p t (k) ~ — — (64) 
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or 



*-<*> = ^ - 1ST {k>>0) (65) 

6 The Sum Rule. 

The crowning result is Levinson's theorem in one dimension: integrating (jfilj) from k = 
to k = oo we find 

# = MO) - yt(oo)] _ 6(0) (66) 

7T 2 

as first given (except that he took only 6(0) = —1) by Barton [TT] . 

We already showed that 6(0) = —1 for a generic potential, that is, when the full 
wavefunction tpk{x) — > for k — > 0. When the potential is critical the zero-energy wave 
function is non-zero, just starting from 

ipk->o = [exp(ikx) + b{k) exp(— ikx)]k=o = 1 (a; < 0) (67) 

and so 6(0) = for u{x) critical. Hence then |£(0)| = 1, but the phase depends on the 
potential. In particular if u(x) is even, the zero-energy resonance is either even or odd and 
therefore £(0) = ±1, ipk=o{x) even/odd. For example, in the Poschl- Teller case given by 



u(x) = — , £ integer (68) 



+ 1) ; 
cosh 2 x' 

the zero-energy resonance has the parity of I. Thus the first member t = 1 of the series 
obtains (see flH| ) 



ipk=o(x) = — tanh x. (69) 

Now we come to the most conspicuous aspect of the one-dimensional scattering, namely 
the one-half factor in the generic case t)66j) . As hinted at by Barton, the reason is related 
to the fact that any attractive potential binds in one-dimension (this theorem seems due 
to R. Peierls ^Hj). Then, the u(x) = potential is critical, that is, increasing it a little bit 
in the attractive side produces a bona fide bound state. Indeed, b(k) = for no potential, 
the earmark for a critical potential, namely transparency at k = 0. We like to call it 
supercritical because it is transparent, i.e. there is no reflection at any energy. 

The phase of the forward amplitude is related to the determinant of the S-matrix. From 
(JSJ) and the phase relation (JTUJ) 



Det S(k) =tt-bb = t 2 -\b\ 2 expi(pr + ¥r) = \t\ 2 exp(2i<^ t ) - e iw \b\ 2 exp(2i<p t ) = exp(2i<p t ) 

(70) 

already noticed by many people, e.g. [0]. 
Therefore, if s(k) := DetS(k), 
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^ N + b{ ^ = ^ll i s dk < 71 > 
very similar to the conventional proof, as the integral can be performed in the complex 
plane [3]. 

The 1/2 contributions for critical potentials, both free u = and interacting, are 
reminiscent of the ^-invariant in the APS index theorem for manifolds with boundary; 
indeed, this can be seen explicitly in the supersymmetric formulation, in which it appears 
as the index of the Dirac operator, giving e.g. fermion numbers 1/2 (see jTOJ and (2D])- 
Moreover, the 1/2 value is characteristic of time- reversal invariant systems, both here and 
in the fractionization case. 

7 Final remarks. 

In principle our method can be applied in arbitrary dimension D: the Schrodinger equation 
for scattering reads 

V 2 V>(r) + k 2 i[)(r) = u(r) i/;(r) (72) 

where r G R D . Again u(f) is real so taking real and imaginary parts and differentiating 
with respect to k, we can get e.g. for the real part $ = 9?{-?/>} 

/ [$V$ - $V$] d D - x o = 2k [ $ 2 d D f (73) 
Jav Jv 

The first term is evaluated asymptotically in terms of 

^k{r) -> exp(ikr cos 6) + r~ (D ~ 1)/2 /(^) exp(i£r) (74) 

as r ^> 0, without supposing u(r) = u(\r\), i.e. non-central potentials are included. The 
second term in (|75j) is related, as before, to the spectral density. The calculation proceeds 
in the same way, except that for D > 2 the case u(r) = is not critical. We refrain to 
reproduce the well known results both for D = 2 and for D > 2. (For D = 2 see ^2], jH], 

i) 

Non-local potentials require a different strategy, because then time reversal T does 
not necessarily hold. Here we want just to show how in one-dimension a non-local real 
potential, which is T invariant, gives, in general, the same result as the local case. The 
Schrodinger equation is now 



ip"(x) + k 2 il)(x) = l u(x,y)tfj(y) dy (75) 

J — oo 

where 

u(x,y) = u(y,x)* (76) 
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in all cases from hermiticity of the hamiltonian H. If moreover T reversal holds, u(x,y) is 
real (hence symmetric) and we get from ()72|) 

/oo 
u(x,y)tp(y) dy (77) 
-oo 

Once again it suffices to take real and imaginary parts and write the wronskian to 
eliminate the term of the potential, so Levinson's theorem seems to hold untouched. Ho- 
wever, for a non-local potential there might be exceptionally bound states embeded in the 
continuum; for the form of the theorem in these cases see |21| . 

As a final comment we want to compare the optical theorem in D dimensions 
with this Levinson's theorem. Both depend on hard data, hence the appearance of the 
forward amplitude is to be expected. Also, they are interference-type formulae, linear on 
t, and represent the same completeness. For the optical theorem it is in coordinate space 
and takes the form of a conserved current, indeed the Noether current associated to the 
classical lagrangian reproducing the time-dependent Schrodinger equation with a global 
phase invariance. For the spectral density the completeness appears in fc-space; the sum 
rule for this case is a kind of global invariant of the problem. 

The generalization of this presentation for the arbitrary D dimensional case is in 
progress 1 . 

Appendix 

Equation is 

I = r *-y- sin((^ r + 2kL) dk for L ->• oo (78) 
Define 2kL = h ! \ then 



/OO POO 
b{k'/2L) sin[(p r {k' /2L) + k') dk' /k' = 6(0) / sin(^) dk' /k' = nb(0) (79) 
-oo J—oo 

as both <f r (0), <^ r (0) and b(0) are regular. 
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